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Nonlinear theory of a "shear-current" effect and mean-field magnetic dynamos 

Igor Rogachevskijj and Nathan Kleeoriijj] 

Department of Mechanical Engineering, 

The Ben-Gurion University of the Negev, 

FOB 653, Beer-Sheva 84105, Israel 

(Dated: February 2, 2008) 

The nonlinear theory of a "shear-current" effect in a nonrotating and nonhelical homogeneous 
turbulence with an imposed mean velocity shear is developed. The "shear-current" effect is associ- 
ated with the WxJ-term in the mean electromotive force and causes the generation of the mean 
magnetic field even in a nonrotating and nonhelical homogeneous turbulence (where W is the mean 
vorticity and J is the mean electric current), ft is found that there is no quenching of the nonlinear 
"shear-current" effect contrary to the quenching of the nonlinear Q-effect, the nonlinear turbulent 
magnetic diffusion, etc. During the nonlinear growth of the mean magnetic field, the "shear-current" 
effect only changes its sign at some value B, of the mean magnetic field. The magnitude B, de- 
termines the level of the saturated mean magnetic field which is less than the equipartition field, 
ft is shown that the background magnetic fluctuations due to the small-scale dynamo enhance the 
"shear-current" effect, and reduce the magnitude B,. When the level of the background magnetic 
fluctuations is larger than i/3 of the kinetic energy of the turbulence, the mean magnetic field can 
be generated due to the "shear-current" effect for an arbitrary exponent of the energy spectrum of 
the velocity fluctuations. 

PACS numbers: 47.65.-(-a; 47.27.-i 
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I. INTRODUCTION 

The magnetic fields of the Sun, solar type stars, galax- 
ies and planets are believed to be generated by a dy- 
namo process (see, e.g., 0,1111111111 111111 111). In 
the framework of the mean-field approach, the large-scale 
magnetic field B is determined by the induction equation 



r)lr\ - - 

^— = Vx (UxB- 
dt ^ 



f:(B)-7?VxB) 



(1) 



where U is the mean velocity, ij is the magnetic diffusion 
due to the electrical conductivity of fiuid. The mean 
electromotive force £(B) = (u x b) is given by 

£,(B) = ay(B)S, - 7y,,(B)(VxB), + (V°ff(B)xB), 
-[^(B)x(VxB)], - K,,k{B){dB),k , (2) 

where u and b are fluctuations of the velocity and mag- 
netic field, respectively, angular brackets denote ensem- 
ble averaging, {dB)ij = C^iBj + \7jBi)/2 is the sym- 
metric part of the gradient tensor of the mean magnetic 
field V^BJ, i.e., V^Bj = {dB),j -1- £„„(VxB)„/2 arid 
Eijk is the Levi-Civita tensor. Here Q:y(B) and ?7ij(B) 
determine the a effect and turbulent magnetic diffu- 
sion, respectively, V'^*^(B) is the effective drift velocity 
of the magnetic field, ^^^(B) describes a contribution to 
the mean electromotive force related with the symmetric 
parts of the gradient tensor of the mean magnetic field. 
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{dB)ij, and arises in an anisotropic turbulence, and the 
^(B)-term determines a nontrivial behavior of the mean 
magnetic field in an anisotropic turbulence. 

The mean magnetic field can be generated in a heli- 
cal rotating turbulence due to the a effect described by 
aij{'B)Bj term in the mean electromotive force. When 
the rotation is a nonuniform, the generation of the mean 
magnetic field is caused by the ail-dynamo. For a rotat- 
ing nonhelical turbulence the 5-term in the mean elec- 
tromotive force describes the fJ x J-effect which causes 
a generation of the mean magnetic field if rotation is a 
nonuniform (see [lU Il3. Ila. Il4l Il5 | ) . where ft is the an- 
gular velocity and J is the mean electric current. 

For a nonrotating and nonhelical turbulence the a ef- 
fect and the fi x J-effect vanish. However, a mean mag- 
netic field can be generated in a nonrotating and nonheli- 
cal turbulence with an imposed mean velocity shear due 
to the " shear-current" effect il9| , described by the ^-term 
in the mean electromotive force. In order to elucidate the 
physics of the "shear-current" effect, we compare the a- 
effect in the arj-dynamo with the ^-term caused by the 
"shear-current" effect. The a-term in the mean electro- 
motive force which is responsible for the generation of 
the mean magnetic field, reads £°' = aBi ex — (J7 • A)B.i 
(see, e.g., d [il|), where A = V(u^)/(u^) determines 
one of the inhomogeneities of the turbulence. The ^-term 
in the electromotive force caused by the "shear-current" 
effect is given by £,f = -(5x(VxB))i ex -(W • V)B, 
(see Eq. (jSJ below, and 16]), where the ^-term is pro- 
portional to the mean vorticity W = VxU. The mean 
vorticity W in the "shear-current" dynamo plays a role 
of a differential rotation and an inhomogencity of the 
mean magnetic field plays a role of the inhomogencity 
of turbulence. During the generation of the mean mag- 
netic field in both cases (in the ail-dynamo and in the 



"shear-current" dynamo), the mean electric current along 
the original mean magnetic field arises. The a-effect 
is related with the hydrodynamic helicity ex ($7 • A) in 
an inhomogeneous turbulence. The deformations of the 
magnetic field lines are caused by upward and downward 
rotating turbulent eddies in the aJl-dynamo. Since the 
turbulence is inhomogeneous (which breaks a symmetry 
between the upward and downward eddies), their total 
effect on the mean magnetic field does not vanish and it 
creates the mean electric current along the original mean 
magnetic field. 

In a turbulent flow with an imposed mean velocity 
shear, the inhomogeneity of the original mean magnetic 
field breaks a symmetry between the influence of upward 
and downward turbulent eddies on the mean magnetic 
field. The deformations of the magnetic field lines in the 
" shear-current" dynamo are caused by upward and down- 
ward turbulent eddies which result in the mean electric 
current along the mean magnetic field and produce the 
magnetic dynamo. 

The " shear-current" effect was studied in [13 in a kine- 
matic approximation. Kinematic dynamo models predict 
a field that grows without limit, and they give no esti- 
mate of the magnitude for the generated magnetic field. 
In order to find the magnitude of the field, the nonlinear 
effects which limit the field growth must be taken into 
account. 

The main goal of this study is to develop a nonlinear 
theory of the "shear-current" effect. We demonstrated 
that the nonlinear "shear-current" effect is very impor- 
tant nonlinearity in a mean-field dynamo. During the 
nonlinear growth of the mean magnetic field, the "shear- 
current" effect changes its sign, but there is no quenching 
of this effect contrary to the quenching of the nonlin- 
ear a-effect, the nonlinear turbulent magnetic diffusion, 
etc. The nonlinear "shear-current" effect determines the 
level of the saturated mean magnetic field. This paper 
is organized as follows. First, we discuss qualitatively 
a mechanism for the "shear-current" effect (Section II). 
In Section III we formulate the governing equations, the 
assumptions and the procedure of the derivation of the 
nonlinear mean electromotive force in a turbulent flow 
with a mean velocity shear. In Section IV we analyze the 
coefficients defining the mean electromotive force for a 
shear-free turbulence and for a sheared turbulence, and 
consider the implications of the obtained results to the 
mean-field magnetic dynamo. The nonlinear saturation 
of the mean magnetic field and astrophysical applications 
of the obtained results are discussed in Section V. 



motive force can be written in the form: 



Ci — Clij Dj -f- O^j}^ ^j^k ) 



(3) 



where Bji = VjBj and we neglected terms ^ 0(V^i?fc). 
Following to [l7J we rewrite Eq. Q for the mean electro- 
motive force in the form of Eq. (0) with 



1. ^ -.off,-. 1 



ay(B) = ^K+OjO , ^r(B)-^£fcjza^j , (4) 

VijK'^) = T(,£ifcp Ujkp + ^jkp Oikp) , Oi — -ryOjji ~ Ojij) , 



(5) 
(6) 



HijkCB) — —-{bijk+bikj), 

where we used an identity Bj^i — {dB)ij + 
£yn(VxB)„/2. Note that the separation of terms in 
Eqs. Q-® is not unique, because a gradient term can al- 
ways be added to the electromotive force. Let us consider 
a homogeneous, nonrotating and nonhelical turbulence. 
Then in the kinematic approximation the tensor a^ van- 
ishes. This implies that a^ = and V^^ = 0. The mean 
electromotive force 5 is a true vector, whereas the mean 
magnetic field B is a pseudo- vector. Thus, the tensor 
bijk is a pseudo-tensor (see Eq. Q). For homogeneous, 
isotropic and nonhelical turbulence the tensor hijk = 
rjrpSijk, where -q^ — uqIq/S is the coefficient of isotropic 
turbulent magnetic diffusion, uq is the characteristic tur- 
bulent velocity in the maximum scale of turbulent mo- 
tions Iq. In a turbulent flow with an imposed mean ve- 
locity shear, the tensor bijk depends on the true tensor 
WjUi. In this case turbulence is anisotropic. The tensor 
Vj Ui can be written as a sum of the symmetric and an- 
tisymmetric parts, i.e., S/jUi — {dU)ij — {l/2)eijkWk, 
where {dU)ij — C^iUj + \7jUi)/2 is the true tensor and 
the mean vorticity W is a pseudo- vector. Now we take 
into account the effect which is linear in \/jUi. Thus, the 
pseudo-tensor bijk in the kinematic approximation has 
the following form 

bijk ~ Vr^ijk + ^0 [^1 £ijm{dU)„ik + D2 £ikm{dU)mj 

+D3 Sjk,n{dU)„u + D4 SijWk + D5 SikWj] , (7) 

where Dk are the unknown coefficients, Sij is the Kro- 
necker tensor, and the term ex 6jkWi vanishes since 
V • B = (see Eq. |j3Jl). Using Eqs. (gJ-jHl we determine 
the turbulent coefficients defining the mean electromo- 
tive force for a homogeneous and nonhelical turbulence 
with a mean velocity shear: 



II. THE "SHEAR-CURRENT" EFFECT 

In order to describe the " shear-current" effect we need 
to determine the mean electromotive force. The general 
form of the mean electromotive force in a turbulent flow 
with a mean velocity shear can be obtained even from 
simple symmetry reasoning. Indeed, the mean electro- 



mj = VT^v-^llvoidU),,, S^l^SoW, (8) 



Kijk = ll [ki 5ij Wk + K2 £ym (9t/)„ifc] , 



(9) 



where 770 = [Di - D2 - 2D3)/i, 5o = (D^ - D^)/2, 
Ki = —{Di + Ds) and K2 = —(Di + 02). The second 
term in the tensor rjij describes an anisotropic part of 
turbulent magnetic diffusion caused by the mean velocity 



shear, while the first term in the tensor rjij is the isotropic 
contribution to turbulent magnetic diffusion. The S term 
for the mean electromotive force describes the "shear- 
current" effect which can cause the mean-field magnetic 
dynamo. Indeed, consider a homogeneous divergence-free 
turbulence with a mean velocity shear, U = (0, Sx, 0) 
and W = (0, 0, S). Let us study a simple case when the 
mean magnetic field is B = {Bx{z),By{z),0). The mean 
magnetic field in the kinematic approximation is deter- 
mined by 



dt 

dBy 

dt 



-SllaoB'^ + V^K, (10) 

5S,+7?^b;', (11) 



where B" = d'^B/dz^ and 



Co = (5o - m 



Hi 
2 



4 



1 



{D2 + Ds + 2D4) . (12) 



In Eq. Hll|l we took into account that the characteristic 
spatial scale Lb of the mean magnetic field variations is 
much larger than the maximum scale of turbulent mo- 
tions Iq. Equation H12(l was obtained in n^ in the kine- 
matic approximation. A solution of Eqs. H1U|) and (|11|1 we 
seek for in the form oc exp{'-ft + iKz z), where the growth 
rate 7 of the magnetic dynamo instability is given by 



J = SIqK^ y/a^ - 7]^ KI 



(13) 



The first term (ex SB^) in RHS of Eq. JTTl) describes the 
shear motions. This effect is similar to the differential 
rotation because Vx(UxB) = SB^Gy The magnetic 
dynamo instability is determined by a coupling between 
the components of the mean magnetic field. In particular, 
the inhomogeneous magnetic field By generates the field 
Bx due to the "shear-current" effect (described by the 
first term in RHS of Eq. (fTn| ). This is similar to the 
a effect. On the other hand, the field Bx generates the 
field By due to the pure shear effect (described by the first 
term in RHS of Eq. (|ll|l ). like the differential rotation in 
afJ-dynamo. It follows from Eqs. (|10|l and (|ll(l that for 
the "shear-current" dynamo, B^/By ^ Iq/Lb <C 1. Note 
that in the aJl-dynamo, the poloidal component of the 
mean magnetic field is much smaller than the toroidal 
field. 

The magnetic dynamo instability due to the "shear- 
current" effect is different from that for afi-dynamo. In- 
deed, the dynamo mechanism due to the " shear-current" 
effect acts even in homogeneous nonhelical turbulence, 
while the alpha effect vanishes for homogeneous turbu- 
lence. 

The "shear-current" effect was studied in 0| in the 
kinematic approximation using two different methods: 
the r-approximation (the Orszag third-order closure 
procedure) and the stochastic calculus (the path in- 
tegral representation of the solution of the induction 
equation, Feynman-Kac formula and Cameron-Martin- 
Girsanov theorem). The 5-term in the electromotive 



force which is responsible for the "shear-current" effect 
was also independently found in [13 in a problem of a 
screw dynamo using the modified second-order correla- 
tion approximation. 



III. THE GOVERNING EQUATIONS AND THE 

PROCEDURE OF THE DERIVATION OF THE 

NONLINEAR EFFECTS 

Now let us develop a nonlinear theory of the "shear- 
current" effect. In order to derive equations for the non- 
linear coefficients defining the mean electromotive force 
in a homogeneous turbulence with a mean velocity shear, 
we will use a mean field approach in which the mag- 
netic and velocity fields are divided into the mean and 
fluctuating parts, where the fluctuating parts have zero 
mean values. The procedure of the derivation of equa- 
tion for the nonlinear mean electromotive force is as fol- 
lows (for details, see Appendix A). We consider the case 
of large hydrodynamic and magnetic Reynolds numbers. 
The momentum equation and the induction equation for 
the turbulent fields are given by 



9u(i,x) 
dt 

dh{t,x) 



_ Vptot + ^[(b . v)B + (B • V)b] 
P MP 

-(U- V)u-(u- V)U + u^ + F, (14) 

(B- V)u-(u- V)B-(U- V)b 

+ (b-V)U + b^, (15) 



where u and b are fiuctuations of velocity and mag- 
netic field, respectively, B is the mean magnetic field, 
U is the mean velocity field, p is the fiuid density, /i 
is the magnetic permeability of the fiuid, F is a ran- 
dom external stirring force, u^ and b^ are the nonlin- 
ear terms which include the molecular dissipative terms, 
Ptot = P + M^^ (B • b) are the fluctuations of total pres- 
sure, p are the fiuctuations of fiuid pressure. The velocity 
u satisfies to the equation: V • u = 0. Hereafter we omit 
H in equations, i.e., we include ^~^''^ in the definition of 
magnetic field. We study the effect of a mean velocity 
shear on the mean electromotive force. 

Using Eqs. (|14() - (|15|l written in a Fourier space we de- 
rive equations for the correlation functions of the veloc- 
ity field fij = {uiUj), the magnetic field hij = (bibj) and 
the cross-helicity gij = (biUj). The equations for these 
correlation functions are given by Eqs. HA1|1 - HA3|1 in Ap- 
pendix A. We split the tensors fij, hij and gij into non- 

helical, fij, and helical, /^ , parts. The helical part of 

the tensor h^- for magnetic fluctuations depends on the 
magnetic helicity, and it is determined by the dynamic 
equation which follows from the ma gnet ic helicity conser- 
vation arguments (see, e.g., [T9II20I mll2l I2ll2^ V The 
characteristic time of evolution of the nonhelical part of 
the magnetic tensor hij is of the order of the turbulent 
correlation time tq = Iq/uq, while the relaxation time of 



,(ff) 



is of the or- 



the helical part of the magnetic tensor h^ 

der of To Rm (see, e.g., 1231), where Rm = IoUq/t] :^ 1 is 

the magnetic Reynolds number, uq is the characteristic 

turbulent velocity in the maximum scale Iq of turbulent 

motions. 

Then we split the nonhelical parts of the correlation 
functions fij , hij and gij into symmetric and antisymmet- 
ric tensors with respect to the wave vector k, e.g., fij = 



As) , Aa) 



p(«) 



f^l' + f^'^ ^here the tensors f^' = [/,,(k) + /.,(-k)]/2 
describes the symmetric part of the tensor and fy^' — 
[/y(k) — /.y(— k)]/2 determines the antisymmetric part 
of the tensor. 

Equations for the second moments contain higher mo- 
ments and a problem of closing the equations for the 
higher moments arises. Various approximate methods 
have been proposed for the solution of problems of 
this type (see, e.g., [23, lla, |23). The simplest pro- 
cedure is the T-approximation, which is widely used in 
the theory of kinetic equations, in passive scalar turbu- 
lence and rna g netohydrod y namic turbulence (see, e.g., 
miii m Ho, m iHl Ha). This procedure allows us 
to express the deviations of the third moments from the 
background turbulence in terms of the corresponding de- 
viations of the second moments, e.g.. 



r - f 



N{0) 



(0)n 



-(/..- i^'OMfc) 



(16) 



where the term ff^ is related with the third moment 
(see Appendix A). The superscript (0) corresponds to 
the background turbulence (with B = 0), and r(fc) is the 
characteristic relaxation time of the statistical moments. 
We applied the r-approximation only for the nonhelical 
part hij of the tensor of magnetic fluctuations. 

In this study we consider an intermediate nonlinearity 
which implies that the mean magnetic field is not enough 
strong in order to affect the correlation time of turbulent 
velocity field. The theory for a very strong mean mag- 
netic field can be corrected after taking into account a 
dependence of the correlation time of the turbulent ve- 
locity field on the mean magnetic field. We assume that 
the characteristic time of variation of the mean magnetic 
field B is substantially larger than the correlation time 
T(fc) for all turbulence scales (which corresponds to the 
mean- field approach). This allows us to get a station- 
ary solution for the equations for the second moments 
fij, hij and gij. For the integration in k-space of these 
second moments we have to specify a model for the back- 
ground turbulence (with B = 0). We use the follow- 
ing model for the background homogeneous and isotropic 
turbulence: fg\k) = (u2)(o) W^(fc) (5,, - fc,,), hfj\k) = 

(b2)(o)T/F(fc)((5y - k,j) and ^^^(k) = 0, where % = 
kikj/P, W{k) = -{df{k)/dk)/8T:P, T{k) = 2Tof(fc), 
f(fc) = (fc/fco)^^', 1 < g < 3 is the exponent of the 
kinetic energy spectrum (e.g., q — 5/3 for Kolmogorov 

f(0) 



Using the derived equations for the second moments 
fij , hij and gij we calculate the mean electromotive force 

£i = / fi(k) dk, where fi(k) = eiTnnffnm(k). For a turbu- 
lence with a mean velocity shear the coefficients defining 
the mean electromotive force are the sum of contribu- 
tions arising from a shear-free turbulence and sheared 
turbulence (see Section IV). 



IV. THE NONLINEAR MEAN-FIELD DYNAMO 

IN A TURBULENCE WITH A MEAN VELOCITY 

SHEAR 

First, let us consider a shear-free nonrotating homo- 
geneous and nonhelical turbulence. Using Eqs. ©-® 
and HA35|I - HA36() we derive equations for the mean elec- 
tromotive force. The coefficients defining the mean elec- 
tromotive force for a shear-free turbulence in a dimen- 
sionless form are given by 



,.(™) 



= a(")(B)%, 



V^« ^ YA{B)+fi{B) 



(B ■ V)B 
B2 



V^j = Va (B) S^j , 



(17) 

(18) 
(19) 



where fi{B) = -(1/2)(1 + e)A)^'{AB) , the functions 
rij^{B) and ~V a{B) are determined by Eqs. H22|) and 

(I24|l respectively, the functions A^, ' (y) are determined 
by Eqs. IJC1|I in Appendix C, the parameter e = 
(b^)(°)/(u^)^°) is the ratio of the magnetic and kinetic 
energies in the background turbulence. The function 
a(™)(B) = x^'^HB) (^("HVSB) is the magnetic part of 
the a-effect, where (j)'^™'\y) — (3/y^)(l — arctany/y) 
is the quenching function of the magnetic part of the 
a-effect (see |3l|, |3J|), and the dimensionless function 
x(^)(B) = (T/3AtpMo)(b • (Vxb)). The function x^^K^) 
is determined by the dynamic equation which follows 
from the magnetic helicity conservation arguments (see, 
e.g., [H IMIH miUlIll). Note that in a homoge- 
neous and nonhelical background turbulence the hydro- 

(u) 

dynamic part, a^,- , of the a effect vanishes. In a tur- 
bulence without a uniform rotation or a mean velocity 
shear, the ^(B)-term and the Kijfc(B)-term in the mean 
electromotive force vanish. 

We adopt here the dimensionless form of the mean dy- 
namo equations; in particular, length is measured in units 
of L, time is measured in units of L"^ /t]^ and B is mea- 
sured in units of the equipartition energy _Boq = y/JIpuo, 
the nonlinear turbulent magnetic diffusion coefficients are 
measured in units of the characteristic value of the turbu- 
lent magnetic diffusivity rj^ = IqUq/S. Note that L ~ Lb, 
where Lb is the characteristic scale of the mean magnetic 
field variations. 

Now we consider a small-scale homogeneous turbulence 



spectrum), fco = l/lo, and tq = lo/uo, J /.^ (k) dk = ^^^j^ ^ ^^^^ ^^^^^-^^ ^^^^^^ U = 5xe, and W = ^e 
((u^)(''^/3)(5ij and //ij-- (k)dk = ((b^)(")/3)(5ij. In cartesian coordinates the mean magnetic field, B 



B(x, z) Gj, + Vx [A{x, z) ey], is determined by the dimen- 
sionless dynamo equations 

^ = (I) S.MB)(W-V)B + a^"^HB)B 

-{VA{B)-V)A + r^,{B)AA, (20) 

— = -5. (W • V) A + V • {r,, {B)VB) , (21) 

where 5* — SL'^/r]^, and W — W /W , the function 
(Jq{B) is determined below, the nonhnear turbulent mag- 
netic diffusion coefficients and the nonlinear drift veloci- 
ties of the mean magnetic field are given by 



i(i 



r,AB) = AY>{AB)+A\'>{AB), 



Wt 



r^,{B) = 4i)(4B) + 3(l-e) 
1 



a^^\ab) 



^a{B) = - 



27r 



^2(165^ 



(2-3e)A«(4B) 



-(1 



e)— A2(16i?^) 

ZTT 



(22) 



(23) 



(24) 



where A^^^ = VB^/B^, the parameter < e < 1, the 

functions Ak{y) and A\. ' (y) are determined by Eqs. IJB6|) 
and (|C1|) in Appendixes B and C. For derivations 
Eqs. (EH)- (EH) we used Eqs. ^ and ^. Note that 
in Eqs. (|22|1 - H24|I we neglected small contributions ^ 
O[{lo/L)'^] caused by the mean velocity shear. The non- 
linear turbulent magnetic diffusion coefficients rj^ and 
Tjg and the nonlinear effective drift velocity Va of mean 
magnetic field for different value of the parameter e are 
shown in FIGS. 1-2. The background magnetic fluctua- 
tions caused by the small-scale dynamo result in increase 
of the nonlinear turbulent magnetic diffusion coefficient 
rjg , and they do not affect the nonlinear turbulent mag- 
netic diffusion coefficient j]^ (see FIG. 1). On the other 
hand, the background magnetic ffuctuations strongly af- 
fect the nonlinear effective drift velocity Va of mean mag- 
netic field. In particular, when e > 1/2, the velocity Va 
is negative (i.e., it is diamagnetic velocity) which causes 
a drift of the magnetic field components B^ and B^ from 
the regions with a high intensity of the mean magnetic 
field B. When < e < 1/2, the effective drift velocity Va 
is paramagnetic velocity for a weak mean magnetic field 
(see FIG. 2). For strong field, B > 5cq/2, the effective 
drift velocity Va is diamagnetic for an arbitrary level of 
the background magnetic fluctuations. 

The asymptotic formulas for the magnetic part of the 
a-effect, the nonlinear turbulent magnetic diffusion co- 
efficients, and the nonlinear drift velocity of the mean 
magnetic field for B ^ -Bcq/4 are given by 



-^(™)(B) = x(^)(B)fl-^ 




3 B/B 



eq 




B/B 



eq 



FIG. 1: The nonlinear turbulent magnetic diffusion coeffi- 
cients T]^ (FIG. la) and rj^ (FIG. lb) for e = (solid) and 
e — 1 (dashed). The function rj^ is independent of the param- 
eter e. The nonlinear turbulent magnetic diffusion coefficients 
are measured in units of the characteristic value of the turbu- 
lent magnetic diffusivity 77^ = Iquq/S. 



VAB) = l-^P\ 77«(B) = l-^(5-4e)/32, 



Va{B) = -(l-2e)/32A(B), 
5 

and for B ^ i?oq/4 they are given by 
VAiB) = -i±^A(^) 



3/3 ' 



3/3 



V 

A 



B <C -Boq/4 is given by 
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FIG. 2: The nonlinear effective drift velocity Va of mean 
magnetic field for e = (solid); e = 0.3 (dashed-dotted); e = 1 
(dashed). The velocity Va is measured in units of r]rp / L. 
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FIG. 3: The dimensionless nonlinear coefficient ao{B) defin- 
ing the "shear-current" effect for different values of the pa- 
rameter e: e = (solid); e = 0.2 (dashed-dotted); e = 1 
(dashed). 



where (3 — ^/SB. 

The nonlinear coefficient (Tq{B) defining the "shear- 
current" effect is determined by Eqs. (|A51|I in Appendix 
A. The nonlinear dependence of the parameter cro(-B) 
is shown in FIG. 3 for different values of the parame- 
ter e. The background magnetic fluctuations caused by 
the small-scale dynamo and described by the parame- 
ter e, increase the parameter ao{B). Note that the pa- 
rameter ao{B) is determined by the contributions from 
the J(B)-term, the 77y(B)-term and the Kijfc(B)-term in 
the mean electromotive force. The asymptotic formula 
for the parameter (7q[B) for a weak mean magnetic field 



ao{B)^-{2^q 



3e) 



(25) 



where q is the exponent of the energy spectrum of the 
background turbulence. In Eq. (|25|l we neglected small 
contribution ^^ 0\(4B /Bpn)^]. Equation (|25|l is in agree- 



ment with that obtained in Il6| where the case of a weak 
mean magnetic field and e = was considered. Thus, 
the mean magnetic field is generated due to the "shear- 
current" effect, when the exponent of the energy spec- 
trum of the velocity fluctuations is 

g < 2 + 3e . 

Note that the parameter q varies in the range 1 < q < 3. 
Therefore, when the level of the background magnetic 
fluctuations caused by the small-scale dynamo is larger 
than 1/3 of the kinetic energy of the velocity fluctua- 
tions, the mean magnetic field can be generated due to 
the "shear-current" effect for an arbitrary exponent q of 
the energy spectrum of the velocity fluctuations. For the 
Kolmogorov turbulence, i.e., when the exponent of the 
energy spectrum of the background turbulence q = 5/3, 
the parameter ao{B) for B ^ -Bcq/4 is given by 



ao{B) = —{l + 9e), 



and for B ^ B^q/A the parameter aQ{B) is 



MB)^-lLil + e) 



(26) 



(27) 



In Eq. H27|) we neglected small contribution ^ 
0{B^q/4:B). It is seeii from Eqs. ^^-^7^ that the non- 
linear coefficient ao{B) defining the "shear-current" ef- 
fect changes its sign at some value of the mean magnetic 
field B = Bf,. For instance, _B» = 0.6i3eq for e = 0, and 
B* = 0.3i?cq for e = 1. The magnitude B* determines 
the level of the saturated mean magnetic field during its 
nonlinear evolution (see Section V). 

Let us determine the threshold for the generation of 
the mean magnetic field due to the "shear-current" ef- 
fect. To this end we introduce the dynamo number in 
the kinematic approximation 



D^ 



loLS 

Vt 



MB^O). 



(28) 



Consider the simple boundary conditions for a layer of 
the thickness 2i in the a;-direction: B{\x\ = l,z) = 
and A{\x\ = l,z) = 0, where x is measured in units L. 
Then Eqs. (HU and ^ yield 



B{t, X, z) 
A{t^x^z) = —Bq 



Boexp{jt) cos{Kxx) cos{Kz z) 



L 



exp(7t) cos{Kxx) sm{Kzz) , 



with the critical dynamo number Da = tt^, where 
ao{B = 0) > 0, the growth rate of the mean mag- 
netic field is 7 = \/D K 



K^ — K^, the wave vector 




0.6 B/B 



eq 



FIG. 4; The normalized nonlinear dynamo number D1f{B) 
for different values of the parameter e: e = (solid); e — 0.2 
(dashed-dotted); e — 1 (dashed). 



K is measured in units of L~^ and the growth rate 7 
is measured in ry^/L^. The mean magnetic field is gen- 
erated when D > D^r- The maximum growth rate of 
the mean magnetic field, 7max ~ -D^/4 — K^, is attained 
at Kz = K,n = vD/2. The critical dynamo number 
determines the critical shear of the mean velocity field 
'S'cr = i''^/i\/o\)){uo/L). The scenario of a nonlinear evo- 
lution of the mean magnetic field is discussed in Section 
V. 



V. DISCUSSION 

In the present paper we developed the nonlinear the- 
ory of the " shear-current" effect in a turbulence with an 
imposed mean velocity shear. The "shear-current" effect 
is associated with the Wx J-term in the mean electromo- 
tive force and causes the generation of the mean magnetic 
field even in a nonrotating and nonhclical homogeneous 
turbulence. The scenario of the mean magnetic field evo- 
lution is as follows. In the kinematic stage, the mean 
magnetic field grows due to the "shear-current" effect 
from a very small seeding magnetic field. During the 
nonlinear growth of the mean magnetic field, the "shear- 
current" effect only changes its sign at some value i?* of 
the mean magnetic field. However, there is no quenching 
of the nonlinear "shear-current" effect contrary to the 
quenching of the nonlinear a-effect, the nonlinear turbu- 
lent magnetic diffusion, etc. The magnitude B^, is less 
than the equipartition field (see below) . The background 
magnetic fluctuations enhance the "shear-current" effect 
and result in a decrease of the magnitude B*. 

The magnitude B* determines the level of the satu- 
rated mean magnetic field. Let us plot the normalized 
nonlinear dynamo number D'^{B) = D'^{B)/ D'^{B — 0) 



which determines the role of the "shear-current" ef- 
fect in the mean magnetic dynamo (see FIG. 4). Here 
D'^{B) = ao{B)/[ri^{B)r]g{B)] is the nonlinear dynamo 
number. At the point B ~ B^ the nonlinear effective 
dynamo number D'^{B) = 0. Depending on the level of 
the background magnetic fluctuations described by the 
parameter e, the saturated mean magnetic field varies 
from 0.3i?cq to 0.6i3eq (see FIG. 4). 

Note that the magnetic part of the a effect caused 
by the magnetic helicity is a purely nonlinear effect. In 
this study we concentrated on the algebraic nonlineari- 
ties (the nonlinear "shear-current" effect, the nonlinear 
turbulent magnetic diffusion, the nonlinear effective drift 
velocity of mean magnetic field) and do not discuss the 
effect of magnetic helicity (the dynamic nonlinearity, see, 
G-g-, [l2,j.2£,.2L^,.23.,.24j) on the nonlinear saturation of 
the mean magnetic field. This is a subject of an ongoing 
separate study. Note that the nonlinear "shear-current" 
effect can affect the flux of magnetic helicity. However, 
this remains an open issue. 

The "shear-current" effect may be important in as- 
trophysical objects like accretion discs where mean ve- 
locity shear comes together with rotation, so that both 
the " shear-current" effect and the a effect might operate. 
Since the " shear-current" effect is not quenched contrary 
to the quenching of the nonlinear a effect, the "shear- 
current" effect might be the only surviving effect, and it 
can explain the dynamics of large-scale magnetic flelds in 
astrophysical bodies with large-scale shearing motions. 



APPENDIX A: THE NONLINEAR MEAN 

ELECTROMOTIVE FORCE IN A TURBULENCE 

WITH A MEAN VELOCITY SHEAR 

We use a mean field approach whereby the velocity, 
pressure and magnetic field are separated into the mean 
and fiuctuating parts, where the fluctuating parts have 
zero mean values. Let us derive equations for the second 
moments. In order to exclude the pressure term from 
the equation of motion (|14|l we calculate Vx(Vxu). 
Then we rewrite the obtained equation and Eq. H15(l in 
a Fourier space. We also apply the two-scale approach, 
e.g., a correlation function 

{ui{x)Uj{y)) = / (ui(ki)uj(k2))exp{i(ki-x 

+k2-y)}dkidk2 
= f Uj (k, K) exp (ik-r + iK-R) dk dK , 



/ fij (k, R) exp (ik-r) dk. 



where hereafter we omitted argument t in the correlation 
functions, /ij(k, R) — L(ui;uj), 



L{a; c) = I (a(k + K/2)c(-k + K/2)) 
X exp (iK-R) dK , 



and R = (x + y)/2, r = x - y, K = ki+k2, k = 
(ki — k2)/2, R and K correspond to the large scales, and 
r and k to the small ones (see, e.g., j^|^). This implies 
that we assumed that there exists a separation of scales, 
i.e., the maximum scale of turbulent motions Iq is much 
smaller then the characteristic scale L of inhomogeneities 
of the mean fields. In particular, this implies that r < 
Iq ^ R- Our final results showed that this assumption 
is indeed valid. We derive equations for the following 
correlation functions: 

/ij(k, R) = L{ui;uj) , hij{k,K) = L{hi\hj) , 
gij{k,R) = L{W,Uj) . 

The equations for these correlation functions are given 
by 



^M^ = ^(k■B)<P^^ - rf 



dt 



dhjj (k) 
dt 



ij '^ ^ijmn\'-')Jmn 



-JM) 



(Al) 



(A2) 



^^^ = ^(k.B)[/,,(k)-;^,,(k)-;.(f)] + 4 



dt 



N 



'■•^ ijmn\^ )9mn + 9ij i 



(A3) 



where hereafter we omitted argument t and R in the cor- 
relation functions and neglected terms ~ O(V^). Here 

$(f )(k) = .9,,(k) - g,,(-k), F.,,{k) = {h{k)u,{^k)) + 

(Mj(k)/'j(-k)), and F(k) = kx(kxF(_k))/fc2p. The ten- 
sors /^,„„(U), E^.^^{\J) and J^,„„(U) are given by 



/"■ (U) 

limn \ J 



d 






^pUq. 



'-^ 'Li'inn \ J 



'^^jq^im^pn ^im^pn^jq ~r OjnOpm^iq 



-rOiinOjn'^q t-.. 



^pUq, 



where Sij is the Kronecker tensor, kij = kikj/k"^. Equa- 
tion (jAl|) - (jA3|) are written in a frame moving with a local 
velocity U of the mean flows. In Eqs. (|AHl and (|A3p we 
neglected small terms which areofthe order of 0(|V^U|). 
The source terms /£ , //i- and If,- which contain the large- 
scale spatial derivatives of the mean magnetic field are 
given by 



tS - 



1 



.iP) 



li^ - -(B.V)$^;'-f[g„(k)(2P,„(fc)-<5™) 

-|-g,,(-k)(2P,„(fc) - (5,„)]B„,, - B„,gfc„$(f^^ 



(A4) 



I^. = -(B.V)<f>(f -b.,(k)J,„-H5,,(-k)5.„]^„., 



R k a)^^-* 
On,qKn^ijg , 



(A5) 






(A6) 



where V = d/OR, <^[p{k) = ^^(k) + g,,(-k), and 
Bij = VjBi, the terms flj, hfj and gfj are determined 
by the third moments appearing due to the nonlinear 

terms, fijq = {l/2)dfij/dkq, and similarly for hijq and 

(p) 
^Ij' . A stirring force in the Navier-Stokes turbulence is 

an external parameter, that determines the background 
turbulence. 

For the derivation of Eqs. IjAip - ljASp we performed sev- 
eral calculations that arc similar to the following, which 
arose, e.g., in computing dgij /dt. The other calculations 
follow similar lines and are not given here. Let us define 
ry(k,R)by 



Fy (k, R) = W {kp + Kp/2)Bp{Q) exp(iK-R) 

x(u,ik + K/2-Q)uj{-k + K/2))dKdCl . (A7) 

Next, we introduce new variables: 

ki = k + K/2-Q, k2 = -k + K/2, (A8) 

k = (ki - k2)/2 - k - Q/2 , K = ki + k2 = K - Q. 

Therefore, 

ry(k,R) = iff,,{k-Q/2,K-Q)ikp + Kp/2)Bp{Q) 
xexp(iK-R)dKdQ . (A9) 

Since |Q| ^ |k| we use the Taylor expansion 
/y(k - Q/2, K - Q) ~ /,,(k, K - Q) 



l df,,{k,K-Q) 
2 dks 

and the following identities: 



Q, + 0(Q2) , (AlO) 



[/,,(k, R)i?p(R)]K = J /y(k, K - Q)Bp{Q) dQ , 

Vp[/,,(k,R)Bp(R)] = J iKp[f,jik,-R)Bp{Ti)]j<, 
X exp (iK-R) dK . (All) 

Therefore, Eqs. (|X9|) - (fnT|) yield 



y,,(k,R) ~ [z(k.B) + (l/2)(B.V)]/,,(k,R) 



1 a^^ 
2^ dK "''' 



(A12) 



We took into account that in Eq. (|A3|I the terms with 
symmetric tensors with respect to the indexes "i" and 
"j" do not contribute to the mean electromotive force 



because £m = £mjigij- In Eqs- (|A9|) - (|A11|I we neglected 
the second and higher derivatives over R. For the deriva- 
tion of Eqs. (|A1(I - (|A3(I we also used the following identity 



iki / /ij(k--Q,K-Q)C/p(Q)exp(iK-R)dKdQ 



1 - 1 - i , - J dfi. 



i fdf,, 



4V dk. 



(V,ViC/p) . 



(A13) 



To derive Eq. (jA13|) we multiply the equation 
V • u = [written in k-space for Wi(ki — Q)] by 
Mj(k2)f/p(Q)exp(iK-R), integrate over K and Q, and 
average over ensemble of velocity fluctuations. Here 
ki = k + K/2 and ka = -k + K/2. This yields 

fi(h + ^K, - QA (u,(k + ^K - Q)uj(-k + ^K)) 
xi7p(Q)exp(iK-R)dKdQ = 0. (A14) 

Now we introduce new variables, ki and k2 determined 
by Eq. HA8|) . This allows us to rewrite Eq. (|A14p in the 
form 



i\h + ^K,~ Q, )/,,(k- iQ,K- Q)c7p(Q) 



2 -^'y^'jv 2 

xexp(iK-R)dKdQ = . 



(A15) 



Since |Q| <C |k| we use the Taylor expansion ljA10|l . and 
we also use the following identities, which are similar to 
Eq. lUTlJ: 

[/,,(k,R)f7p(R)]K = y/^,(k,K-Q)C7p(Q)dQ , 

Vp[/,,(k,R)t/p(R)] = y"^Xp[/„■(k,R)f7p(R)]K 
X exp (iK-R) dK . (A16) 

Therefore, Eq. (|XT5l) yields Eq. (|JT3|) . 

Now we split all tensors into nonhelical, fij, and heli- 



time, and no slightly perturbed steady state can be ap- 
proached |27| . In the r-approximation, the relaxation 
time for small departures from equilibrium is determined 
by the random motions in the equilibrium state, but 
not by the departure from equilibrium ^27] . As follows 
from the analysis by |27l | the r-approximation describes 
the relaxation to equilibrium state (the background tur- 
bulence) much more accurately than the EDQNM ap- 
proach. 



1. Shear-free homogeneous turbulence 



Consider a turbulence without a mean velocity shear, 
i.e., we omit tensors /5,„„(U), £;g„„(U) and J^„(U) 
in Eqs. HA1|I - HA3|I . First we solve Eqs. HAip - ljAsp ne- 
glecting the sources ifj , I^j , Ifj with the large-scale spa- 
tial derivatives. Then we will take into account the terms 
with the large-scale spatial derivatives by perturbations. 
We start with Eqs. ljAip - ljA3|l written for nonhelical parts 
of the tensors, and then consider Eqs. (|A1|) - (|A3|) for he- 
lical parts of the tensors. 

Thus, we subtract Eqs. ljAip - ljA3p written for back- 
ground turbulence (for B = 0) from those for B 7^ 0. 
Then we use the r approximation and neglect the terms 
with the large-scale spatial derivatives. Next, we assume 
that rjk'^ <C t~^ and vk^ -^ t~^ for the inertial range 
of turbulent flow, and we also assume that the charac- 
teristic time of variation of the mean magnetic field B 
is substantially larger than the correlation time T(fc) for 
all turbulence scales. Thus, we arrive to the following 
steady-state solution of the obtained equations: 



3'u(k) 



/(°)(k)+*r(k.B)|.f)(k), 
*r(k-B)[4-(k)-/i,,(k)], 



(A17) 
(A18) 
(A19) 



where fij , hij and gij are solutions without the sources 

Now we split all correlation functions into symmetric 

and antisymmetric parts with respect to the wave number 

cal, f^f\ parts. Note that the helical part of the tensor of k, e.g., f,j ^ f ^f + f^f , where f^f = [f,j{k) + f,j{-k)]/2 

is the symmetric part and f>j = [fij(k) — fij{—k)]/2 is 



magnetic fluctuations h^- depends on the magnetic he- 
licity and it is not determined by Eq. ljA2p . The equation 
for the helical part of the tensor of magnetic fluctuations 
hi, follows from the magnetic helicity conservation ar- 
guments (see, e.g., [H ^ IH HI IH HI ) . 

In this study we use the r approximation [see Eq. 
(|16|l ] . The r-approximation is in general similar to Eddy 
Damped Quasi Normal Markovian (EDQNM) approxi- 
mation. However some principle difference exists between 
these two approaches (see [II i^). The EDQNM clo- 
sures do not relax to equilibrium, and this procedure does 
not describe properly the motions in the equilibrium state 
in contrast to the r-approximation. Within the EDQNM 
theory, there is no dynamically determined relaxation 



the antisymmetric part, and similarly for other tensors. 
Thus, Eqs. (|IT7|) - (|Il9)l yield 



/S^(k) ^ 


\(^ I I'Af'^^Uk) 1 I'-h'^^^ (k)] 


1 + 2i1j ^^' ■'' ' 




(A20) 


/^!;^(k) ^ 


' Y^[^/i^^(k) + (l + ^K°'(k)], 




(A21) 


3?(k) ^ 


^ T+J?[i?^W C(k)]. (A22) 


re V(k) 


= 2(rk-B)^. The correlation fimctions 
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c(a)i 



■M)i 



As) 



f-, ' (k) , hlj (k) and gl- ' (k) vanish if we neglect the large- 
scale spatial derivatives, i.e., they are proportional to 
the first-order spatial derivatives. Equations l|A20l) and 
(|X2T1) yield 



/^;^(k) 



■Mf 



f(0) 



C(k)«C(k) + C(k) 



(0), 



(A23) 



which is in agreement with that a uniform mean mag- 
netic field performs no work on the turbulence. A uni- 
form mean magnetic field can only redistribute the en- 
ergy between hydrodynamic fluctuations and magnetic 
fluctuation. A change of the total energy of fluctuations 
is caused by a nonuniform mean magnetic fleld. 

Next, we take into account the large-scale spatial 
derivatives in Eqs. (jAl|) - (jA3|l by perturbations. Their 
effect determines the following steady-state equations for 
the second moments /^ , hij and ijij : 



M'i^ 



zr(k.B)$(f-)(k) 



?(a) 
lij 

h^\k) = -*r(k.B)$(f-)(k) 
^l-^k) 



tI: 



f 



tL". 



:>lj 



^r(k.B)(i^(;)(k)-/i(;)(k))-K<. 



(A24) 
(A25) 

(A26) 



where l>,^f '"^ = [$»7'(k) + $,^f ^(-k)]/2. The solution 



of Eqs. (IA24I 



$ 



(M,: 



'(k) 



(M), 

26)) yield 

J* 2J / J 






1 



31 



■ZT(k.B)(Z 



(A27) 



Substituting Eq. HA27P into Eqs. (|A24p - HA26ll we obtain 
the flnal expressions in k-space for the nonhelical parts 



c(a)l 



-M 



of the tensors /^"^(k), h\^>{k), gl->{k) and ^\J' '(k). In 
particular. 



.{M,s), 



El m m 113). The solution of Eqs 
yield 



m and JA30 



$ 



(MM) 



(k) 



2ir(k^ (H) 

' 1 + ij '^ 



(A31) 



Since h], is of the order of 0(V) we do not need to take 
into account the source terms with the large-scale spatial 
derivatives pa |. 

Now we calculate the mean electromotive force £i (r — 

0) = (l/2)£„™ J[<^';i'r^"\k) + $L*^'^^(k)] dk. Thus, 



l + 2ij 
»^(k-B) (g 
1 + V "" 



[/f„„-fzT(k.B)(// 

dk. 



mn ran J 



(A32) 



For the integration in k-space of the mean electromo- 
tive force we have to specify a model for the background 
turbulence (with B = 0), see Section III. After the inte- 
gration in k-space we obtain Si = aijBj + bijkBj^ki where 



-le,, 



(H) 



Tfhj flmn 

i + V 



+ 2(3-(l + /3'' 



dk = x^^HB) 



{P)]P^A(3) 



a(™) 



t'^jk = t;Vt 



(A33) 



(1 + £)£,,„ i^6k„rKJ,]}{V2(3) - Kill{V2f3)j 



+2e,nk (1 - e)^i{K,n} - K]l>{V2f3) 



(A34) 



mn \. ) (1+27/^)2 



{l+e){l + 2llj)(5njSmk 



OijijOnk \ '^nk^rnj '^mk^nj ) ■^t,^ ~r '^W jx^^iij^mk 



f£riij ^nk ) 



B 



J,k ■ 



ns), 



;(«)/ 



(A28) 



;^('')f 



The correlation functions f^j'{k), hlj(k) and g,-, (k) 
are of the order of ^ 0{\/'^), i.e., they are proportional to 
the second-order spatial derivatives. Thus fij + fij is the 
nonhelical part of the correlation functions for a shear- 
free turbulence, and similarly for other second moments. 

Now we solve Eqs. (jAl|l - (jA3|l for helical parts of the 
tensors for a shear-free turbulence using the same ap- 
proach which we used in this section. The steady-state 
solution of Eqs. (|A1|) and ljA3p for the helical parts of the 
tensors reads: 



fW, 



fi;>{k) « *r(k.B)$,,^ 



(M,H) 



(k), 



(A29) 



(H) 

9ij 



(k) 



.iM.H) 



ZT(k.B)[/if(k) 



,{H) 



/^(f(k)], (A30) 

AH) 



where ^.^'"'(k) = g^^f ^(k)-^}^ ^(-k). The tensor /i^ 
is determined by the dynamic equation (see, e.g., [13, HO, 



^,{X} = 3X«(V2/3)-:fX(2/32), 



3 

2^ 



^S^W = iz I K,,{a{f))fdf 



3/3* 



K,,{X^ 



X^ 



■dX 



and all calculations are made for q = 5/3, X^ — 
/32(fc//co)2/3 == a = [PuokT{k)/2]^, the function K,^ 
is determined by Eq. (|B1|I in Appendix B, and e — 
(b2)(o)/(u2)(o), (3 = 4B/(uoV2^), P,,(/3) = S., - A,, 
A, =. B,Bj/B\ 0('")(/3) = (3//32)(l - arctan(/3)//3), 
x'^-'(B) = (T/3/xp)(b • (Vxb)) is related with current 
helicity. Since a part of the mean electromotive force is 
determined by the function aij{'B)Bj and Pij{(3)Bj = 0, 
we can drop the term ex Pij (/?) in Eq. HA33(I . Thus, the 
equations for a,jj and bijk are given by 



a(™nB)<5„-, 



6yfc = Vt: 



bl Sijk + &2 Sijn Pnk + ^3 ^ink A 



'ry 



(A35) 



(A36) 
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where ^^(B) = x^'^HB) (/-("H/^), and 

61 = 4'^(V2/3)+A(^)(V2/3), 

b2 = -i(l + e)A«(V2/3), 

63 = (l-e)*i{A2}-4^'(V2/3), 

the functions Ak{y) and A^ (y) are determined by 
Eqs. (jB6p and (|Cip in Appendixes B and C. Equa- 
tions l|J35)l and (|X36|) yield Eqs. (fT7|) - (P|l . 



2. Turbulence with a mean velocity shear 

Now we study the effect of a mean velocity shear on 
the mean electromotive force. We take into account the 
tensors /5„„(U), £^5„„(U) and J5„„(U) in Eqs. (EU- 
(|A3(I . and we neglect terms ~ O(V^). The steady-state 
solution of Eqs. IJA1|) - (JA3|) for the nonhelical parts of the 
tensors for a sheared turbulence reads: 



A^- 



(U)3n 



ijmn 
ijmn 



r{-^(k.B)a>(f)-|-/H 



V 



T{z(k.B)[/,,(k: 



%(k)] 



(A37) 
(A38) 



27 J 



(A39) 



where 



A^/ 



ijmn 
rh 



(U) 



A^" (U) 

ijmn \^ J 



(^imOjn ^ '''-'-ijmn J 
(^imOjn ~ T-t^ijmn 1 
O'imOjn ~ T J ijmn i 



and we use the following notations: the total correlation 
function is /^ = /y -H /!^ . Here /^ = /„■ -h /„ is the 
correlation functions for a shear-free turbulence, and the 
correlation functions f[, determines the effect of a mean 



velocity shear. The similar notations are for other cor- 
relation functions. Now we solve Eqs. (|A37|1 - HA39|I by 
iterations. This yields 



h^^ik) = r{i?5„„ C„ - ^(k•B)ci>, 

9ij['^) — '''{•hjmn9mn 



(M,(t) 



rif:'^) 



}, 



(M,(t) 



4f"^} 



z(k.B)[/--/.^]+4f^^)} 



(A40) 

(A41) 

'}, 
(A42) 



-AM, a) 



where $,■,- ' (k) = fffw(k) — g^j(— k), the source terms 



^y - 4(55)' ^j'''"' = 4i9W and i,^ 



r(/.'^) _ 



{a,<^) - 



Ifjif°j,hfj) are determined by Eqs. (|A4|) - (|A6|) . where /^ , 
/lij, gij are replaced b y f[^, hf,, gfj, respectively. The 
solution of Eqs. 

metric part <& 



(A/,(T,; 



42|l yield equation for the sym- 



$ 



{M,a,s} 



i^)-J 



2V' 
-/if''^)+zr(k.B)[(/^., 



of the tensor 



ijnin ^ jimn J^nin 1 -^^j 






jimnJ J " 



r(/,'^) 



,(/i,^) 



rCi,'^; 



]} 



(A43) 



where we took into account that E",^^ is a symmetric 
tensor in indexes i and j. Thus, the effect of a mean 
velocity shear on the mean electromotive force, £[{r = 
0) = (l/2)£i„,„ / ^Ln dk, is determined by 






1 + 2V''^"'"^'^ 

' ran rnn J ' ran i ^ ' 



X'J mnpq9pq ' ''''''{'^''^JV^mnpqlpq 

(A44) 



Now we use the following identities: 



£impKjkqm^pUq - -[Ci(25j^,, + 25'^^,, + S,^.f, + S'j^ ,, ) + 6*3(25'^^.^, -I- S'j^^,)] , Ei-jmKkmpq^pUq - 2CiS,^.^ 



.(6)^ 



Sivaqi^mj^ kUq = ^ijmi^mq^ kU q = — Ai(ZO^^^, — Jj^jf.) j ^ikqJ^mm^ jU q = —(0^1 -f A2 j (20j^-j, — Jy-j.) , 
SimkKmqW jUq = —£imqKmk^ ]Uq = —-Al{2S^^j^ — S'^^-j.) , Ei-jpRkq^/ pUq — -Al{2S-^f. + S^j^.) , 



1 



ey.i^mmVfcC/, = -(3Ai + A2)(25,;;^ - S'l^D , (£„,j4p + £^,p4g)^mmVpC/, = 2(3^1 + A2)Sl^l , (A45) 



where 
^ijk — £ijp{dU)pk , S^.f. — eikp{dU)pj , 



SSI = ejkp{dU)p^ , Sg = WkS^, , ^,^., 

SSl = e.kp(3jq{dU)pq , 5g = Wk(3^, . 



^S = WjS^k 



12 



After the integration in Eq. (| A44|) . we obtain 
^1 = KikBj,k , 



(A46) 



the coefficient D3 = 0, and the other coefficients calcu- 
lated for q = 5/3 are given by 



where the tensor 6^^, is given by 



b'^ — /2 

"ijk — '0 



r 7 



(A47) 



D^ = 



A^^ - 34'^ - 18Cf ^ + ef^f ^ + 4'^ + ^^f ^) + *i|^i + 2^2 + yCi - e(2Ai + A2 + 6C1) 
+*2|-Ai + \Cx + e^Ai - yCi") I - (1 - e)*3{Ci} - *o{2Ai - 3Ci} 



i.. ^i 



-(Af ) + 4^) + 4Cf' ) + 6( -Af + A'^) + |cf M + * J -Ai + A2 + ^Ci - 2e{A, + 19C,) 



-*2Mi-yGi-e(^i-yC'i 
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(1 - e) ( *3{-2Ai + 7Ci} - y *4{Ci} + 8*5{Ci} 



+ *0<| 2Al - yCl 



i?4 - - 



i^5 



3Af ) + 4') - Hcf ) + e(3Af ) - 4') - |cf )) - ^A, +A,- ^C, -2e[A,+A, 



\Ci 



l-e){ *2{Ai + Ci} - *3{C7i} + *o{C7i} 



4^)+Af 
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<(2) 



l(2) .(2) 



26 



3 ^ V ^ 2 3 



Cf )") - *i|ai - A2 - ^Ci - 2t(A^ - A2 



Ci]\ + {l-e)[ *2{Ai + Ci} - *3{Ci} + *o{Ci} 



-Dfi 



4^)-4Cf-.(Af-|cf 



70 
*!<! -3^2 + yC73 + 2e(A2 - I9C3; 



-(71-79e)*2{C3} 



-(1 - e)( *3{^2 - 7C3} + y*4{C3} - S^sia^}") + *o|^2 - yC3 



^ 6 



4^)-Hcf+/34^)-|c7f 



*3{^2 + C3} + *o{A2 + Ca} 



*i<j A2 + -(1 + e)C3 \ + {l-e)[ *2{2A2 + C3} 



(A48) 



The functions Ak{y) and Ck{y) are determined by ^ rT^i _ ^X^'^''(\/2B) 



(2), 



Eqs. I|B6|I in Appendix B, and the functions A^. ' [y) and 

Cf. (y) are determined by Eqs. IjDljl in Appendix D. The 
functions ^fe{X} are given by 



*o{X} = -i(l + e)x(2)(0) + (2-e)x(2)(y2/3) 



An 



(1 - e)X{2p') , 



*i{X} = -3X(2)(V2/3) + — X(2/32), 



47r 



*3{X} = -6X(2)(V2/3) + 



47r 



*4{^} - 4X(2)(V2/3)--^ 

in 



2X(y)+yX'(y) 

AX{y)+yX'{y) 
SX{y)+AyX'{y) 



y=20^ 



y=20^ 



y'x"{y) 



y=20^ 



^,{X} = -ix(2)(V2/3) + ^ 



2X{y)+yX'{y) 
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+V^"(y) 



(A49) 



y=2/32 



In Eqs. I|A45|I - I|A49|I we took into account that for the 
"shear-current" dynamo, B^/By ^ Iq/Lb ^ 1, where 
Lb is the characteristic scale of the mean magnetic field 
variations. The nonlinear coefficient defining the "shear- 
current" effect is determined by 



MB) 



1 



{D2 + 2Di + D6 + 2Dr) 



(A50) 



Equation HA50I) yields 

ao{B) ^ MAi + A2} + h{Ci + Ca} , (A51) 
where 



l{^} 



(l + 6)X(2)(V2/3) + [*o-(l-e)(*i 



'*2 + *3)]{^} 



(A52) 



h{X} = 9 (3e-13)A(2) (72/3) + [4*2 -4*0 
-18*1 + (1 - e)(55*i - 38*2 + 8*3 
-8*4 + 12*5)]{X} . (A53) 



The nonlinear dependence of the parameter ctq^B) de- 
termined by Eq. ljA5ip . is shown in FIG. 3 for different 
values of the parameter e. The asymptotic formula for 
the parameter (To(-B) for B <C -Bcq/4 and B ^ -Bcq/4 are 
given by Eqs. lf^ - lf?7|) . For the derivation of Eq. (|A5ip 
we used identities (|D2ll in Appendix D. 



APPENDIX B: THE IDENTITIES USED FOR 
THE INTEGRATION IN k-SPACE 

To integrate over the angles in k-space we used the 
following identities: 



K^J = 



^^^■™^ ;ded^^A,s,j+A2p.. 



K, 



1 + a cos^ 9 

^ijmn Sin [ 



tj , 



d9 dip — Ci{6ijSjnn + SiniOjn + OinOjm) + C2pijmn + Cs{dij f3mn + Oimf^jn 



limn — / 1 , 9/1 

■^ J 1 + a COS^ 6 

i^inPjm ~r (JjmPin ~r (JjnPirn ~r OrnnPij ) i 



^ijmn\^} 



'^ijmn Sm ' 



(1 + acos' 

^ijnin Sin ( 



2a-)2 



d9 dip 



d 



^ijmn Sm t 



db J b + a cos^ 



d9 dip 



fc=i 



(Bl) 



(B2) 



Kijmn{a) + a — Ki-imn{a) , (B3) 



(1 -f acos^y)-^ 2\obJ (b + acos-' O)-' J ^_, 2 da 



b=l 



f^ijmnSlllu ^fl ^„ I I "'ijmn Smc 



■ dO dp = 



{l + acos^er-^-^- 3\dbJ (6 + „eos2 0)3'^^^^j,^,=^'^^-(«) + ^-^'^-(")' (^^) 



a d 
3 da 



where a = [/3uoA:T(fc)/2]2, /3,; = /3i//3, /3y = /3i/3j, and 



Al 


2tt ' 
a 


^^^^^arctan(V^) 1 


5 




A. 


27r 


r arctan(V^) 1 






a 


L \/a J 




Ci 


■K 

2a2 


^ ^ _i^2^i"ctan(Y^) 5a _| 
(•^ + Ij 1— -'^ 


C2 


= A2 - 7Ai -1- 35(71 , 


C3 


^ M- 


5Ci . 







In the case of a <C 1 these functions are given by 

T / N 47r/ 1 \ - , , 87r 

AM ^ —\\--ay A^^^ — a, 



(B6) 



^^ ' 105 



In the case of a » 1 these functions are given by 



CM 
Csia) 



TT^ 4:TT T , N TT^ 87r 

—^ , ^2(0) ^H , 

y/a a ya a 

7r2 47r - , , 3tt^ 32tt 

, (--2(0) 



Ay/a 3a 
7r2 87r 



4y^ 3a 



Ay/a 3a 
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APPENDIX C: THE FUNCTIONS aL^^(/?) AND 



The functions An ' (/3) are defined as 



4'H/3) = 



3/34 



An{X^ 

X^ 



■dX , 



and similarly for Cn (/3), where 

X2 == /32(fc/fco)2/3 = a = [/3uofcT(fc)/2]2 , 

and we took into account that the inertial range of the 
turbulence exists in the scales: Id 1^ i" < lo- Here the 
maximum scale of the turbulence Iq <C Lb, and Id — 
lo/Re^''^ is the viscous scale of turbulence, and Lb is the 
characteristic scale of variations of the nonuniform mean 
magnetic field. For very large Reynolds numbers kd = 
l~^ is very large and the turbulent hydrodynamic and 
magnetic energies are very small in the viscous dissipative 
range of the turbulence Q < r < Id- Thus we integrated in 
An over k from feg 
identity 

nl 



L to 00. We also used the following 



2tt 



An{a{f))fdf^—A[lHl3) 



.(1) 



(1)/ 



and similarly for C„ (/3) . The functions Aa (/3) and 
Ca (/3) are given by 



A«(/3) 



A^ 



(1)^ 



2 (/?) = - 



&^\P) 



3 

10 



arctan /3 

~~P 
arctan /3 

. T~ 

arctan /3 



1 



/3 



5 

15 

"^7/32 
10 
^7/32 



rSm - ^ 



7 
5 



il«« 



235 



5 
9^ 



189/32 

A«(/3)-7A«(/3)+35Cf'(/3) 

A«(/3)-5C(^)(/3), 

where L(/3) = 1 - 2/32 + 2/3'* ln(l + /3-2). por /3 < 1 these 
functions are given by 






cf\p) 



1 



:P\ A^i\P) 



:P' 



5V 7'^ 



cf)(/3) 



■S^^^^^' 



35 



P' 



and for /3 3> 1 they are given by 



1(1) 



AY'iP) 
C«(/3) 



5^ " ^ ' ^2 (/3) - - 


37r 4 


20/3 3/32 ' ^2 (P) 


97r 
20/3 ' 


37r 4 




20/3 ' 3/32 ■ 





Here we used that for /3 <C 1 the function i(/3) ^ 1 — 
2/32 - 4/3"* ln/3, and for /3 > 1 the function L(/3) - 2/3/32. 
We also use the identity: 



Hijmn{a{T)) ■ 



df = 27ri^i;i„(/3) 



APPENDIX D: THE FUNCTIONS Ai^'(/3) AND 
Ci'>(/3) 

The functions An (/?) are defined as 



^(2 



and similarly for C„ (/3). We used the following identity 



An{a{f))fUf=^A(^HP), 



and similarly for Cn {P). The functions Aq (/3) and 
Ci (/3) are given by 

4')(/3) = ^(1;-1;0), 

Af\p) = ^(-1;3;0), 

cf)(/3) = (1/4)^(1; -2; 1), 

&i\p) = (1/4)^(3; -30; 35), 

Cf(/3) = (l/4)^(-l;6;-5), (Dl) 

where 

F{a- a; 7) = 7r[a j'"^ (/3) + aj!^'^ (/3) + 7 jf^ (/3)] , 



(CI) jf^(/3) = ^(1 + 



arctan /3 3/32 
^ 2" 



L{P) , 



(2),^^ 7 (2) 



^r(/3) - gjr(/3)+i(/3): 



9 

IT 



jf(/3) = -[j^2\P)-j,m g^^2 



^ . . .. 4 

?2' 



L(/3) 



2 /. arctan/3 2^ 



1 



37r/32 \^^ /3 

For /? <C 1 the functions Ji ^ (/3) are given by 

r(2) 

1 f. 15 






'''^^^-i-A'-T/' 



and for /3 3> 1 they are given by 



jI^\p) 






3 

7/3 
3 


3 

47r/32 


47r/32 
1 


7 



47r/32 



For /3 <C 1 the functions ^L ' (/?) and Ci (/?) are given 

by 



cf\p) 
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14' 



/3^), Cf (/3)^0(/34), 



^(2), 



c^r^(/3) - -^/3' 



and for /3 3> 1 they are given by 

cf(/3) 



cf)(/3) 



Stt 3 

7^~ 2^ ' 
37r 1 

28^" 2/32 ' 
37r 1 

"28/3 ^/32 



4'^(/3) 



37r 3 

"7^5^/32 ' 
97r 4 
28^~^ ' 
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We also used the foUowing identities: 






K^jLiV2P) - H.^ 



SL{V2(3) 



= K^IjV2P)^2hI'IjV2P) 



-G 



(2) 



.(^/2/3) , 



tjmn V 

^SJ^^/3)-2GgL(V2/?) 



■' ijmn V 



<ijV2/3)-3i/i-L(V2;3) 



ijmn \ 



3G\'ljV2P) - Q^'IJV2P) 



ijmn \ 



(D2) 



where 



H. 



(2) 

ij 17171 



(^/2/3) 



4K 



Gl-i„(V2/3) 



(2) 

ijmn 

(2 



(V2/3) 



rj ijmn 



(V2/3) 



gl^L(^^/3) 



3 

2^ 
3 

47r' 
1 

2^' 



i^,- ^-«2 



ijmn\/^P } 



Hi 



.(2/3' 



2G(^L(^^/9)-:iG.,™„(2/32). 
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